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I. Introduction 


This report was prepared In support of the Yield Model Development (YMD) 
Project of the Agriculture and Resources Inventory Surveys though Aerospace 
Remote Sensing (AgRISTARS) program. 

In the YMD project, crop yield models will be used for large area yield fore- 
casting In the U.S. and In foreigi areas both for midseason forecasts and for 
end-of-season estimates of crop yields. Numerous crop growth and yield models 
have been developed for a wide range of purposes in recent years. These models 
range in complexity from the most sophisticated simulators of plant growth, pri- 
marily intended for research into plant physiological interactions, to multiple 
regression models using only a few monthly weather variables to forecast 
regional crop yields. 

In this report we will review plant-process models of a moderate level of 
complexity as candidates for application to large area yield forecasting. 
Identification and selection of such models presupposes some accepted standards 
for model evaluation. However, because of wide differences in philosophy among 
researchers working in plant-process yield modeling, complete agreement on such 
standards has not been reached. Therefore, we will propose a working set of 
standards for identifying and selecting plant-process yield models with poten- 
tial for large area yield forecasting. We will then review some available 
models in terms of these criteria. 

Other criteria were established by Wilson et al., (1980) for selecting sta- 
tistically based crop yield models. The more quantitative of those criteria 
will also be used, when it is appropriate, for evaluation and selection of pro- 
cess based crop yield models. 

II. Plant-Process Yield Modeling 

Numerous plant processes have some degree of importance in determining crop 
yield variability under some range of environmental conditions. For large area 
yield forecasting we want those processes which can capture extreme variations 
which can occur over large areas. 

The processes most basic to yield are photosynthesis and dry matter 
partitioning, as these together directly account for grain formation. 

Phenology, leaf area expansion, florae development, and leaf senescence must 
be accounted for to accurately model photosynthesis and dry matter partitioning 
throughout a crop growth cycle. The number of grains per unit area, which comes 
from floral development, is needed to properly model the partitioning of dry 
matter to grain. Plant water status directly affects all the above processes 
and may be estimated from a soil water budget involving evapotranspiration and 
available moisture capacity of the soil. 

Since these processes are highly Interrelated, one may select a variety of 
different combinations of processes in modeling crop yields. Fbr example, while 
one researcher might directly model photosynthesis, another might try to account 
for its effects by modeling tiller development, leaf area expansion, floral 
development, and water stress. 
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Generally plant-process yield models have been developed to predict yield 
at the level of an average plant In a specified field. Thus the Input data 
required by these models Include plant parameters specific to the variety or 
hybrid planted In some field and soils parameters describing the soli In that 
field. 


Fbr large -area yield forecasting applications, requiring such varietal - 
and field-specific Input data is not practical. This is because within the size 
of areas for which yield predictions may be made (approximately 2,500,000 
hectares) each of several varieties of a major crop will be present on several 
soil types. Therefore plant and soil parameters which represent all or a 
substantial portion of the varieties or hybrids and of the soil types present 
in an area are required. 

III. Model Selection Criteria 

Models will be identified and evaluated. Some will be selected for imme- 
diate testing while others will be modified and eventually tested. Other models 
may supply process submodels to a new model. Models will be selected for 
testing on the basis of how well they meet the needs of the Yield Model 
Development project in the Judgement of project personnel. To aid in that 
Judgement some working standards are proposed. These standards will be used to 
identify plant-process yield models to be acquired for our purposes: for com- 
parative testing; for further development and testing; for incorporation into 
large area yield estimation routines; or to supply subroutines to other models. 

Other criteria may be added as additional insight is gained in this study. 
If any suitable candidate models are not included in this report, it is by over- 
sight and the author requests any Information available on such models. 

Major considerations should Include (1) theoretical accuracy, (2) 
completeness, (3) simplicity, (4) sophistication , (5) structure, (6) validation, 
and (7) timeliness. These will be described in detail in the following 
paragraphs . Additional considerations include public availability of computer 
code, geographic area of model applicability, probable accuracy of model 
estimates, applicability to large area yield forecasting, absence of subjective 
factors in a model, ability to acquire all necessary inputs for a model, and the 
cost of running a model. 

Theoretical Accuracy . Theoretical accuracy refers to how well the model 
agrees with current plant physiological theory. Each component of a model 
should be in agreement with our best current theoretical understanding of how 
that plant process really works. Interactions between model components should 
reflect or at least, not contradict what Is known about how those processes 
interact in real plants. Theoretical accuracy may be assessed by posing the 
following questions: (a) to what extent is the final yield (or an important 
measurable intermediate variable) in the model affected by the full range of 
variation of each submodel or input variable, (b) how much variability in pre- 
dicted yield would be lost by compressing or simplifying any given submodel (or 
variable) into a simple equation (or variable), and (c) how frequent and impor- 
tant are the conditions under which any given submodel will respond differently 
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from the way a real plant responds . While a high degree of theoretical accuracy 
should Increase the probability of a model performing well outside the set of 
conditions where it has been developed and tested, such a model will not 
necessarily give highly accurate yield estimates. This can only be determined 
by a testing program. 

Completeness . Completeness of a model refers to how well that model is 
designed to meet its purpose or objectives considering the limitations under 
which it is to be used. This means that a model should include all those pro- 
cesses and enviroi .mental conditions which are critical to yield formation 
throughout a large area. A model may also include processes which do not 
us uall y directly affect yield but which affect critical processes under certain 
conditions and for which modeling theory is well developed. A model need not 
include processes or conditions for which required parameters or data are not 
readily available, if they affect yield for only a very small part of a large 
area. 

Simplicity . Simplicity is a desirable characteristic of a model. 

Simplicity In model form and use of input data, and the availability of model 
input data are often assoc. iated with cost. Another very Important advantage of 
model simplicity is the ability of the user to understand the concept, capabili- 
ties and limitations of the model. A thorough understanding allows the user to 
evaluate the model in the light of other Information and to modify the model if 
necessary. A simple model will generally require less training and experience 
on the part of the user. 

Sophistication . The sophistication of a model relates to the complexity of 
interactions between model processes and the degree of detail with which each 
process is simulated. A high degree of sophistication can make a model sen- 
sitive to the effects of many combinations of conditions. However sophistica- 
tion can also reduce a model’s simplicity and make the model difficult to 
understand and expensive to run. A model may also be unsuitable for large area 
yield forecasting if input data are required that are not readily available. In 
general a modeler nust find a balance between simplicity and sophistication that 
is appropriate for a particular application. 

Structure. Structure refers to how a model has been Implemented In computer 
code. A properly structured model will be relatively easy to understand and 
needed changes (for evaluation, further development, or operational 
implementation) can be made quickly and accurately only if the model program has 
modularity (black box structure). An Improperly structured model is practically 
impenetrable to everyone but Its originators and proper documentation is all but 
Imixossible. 

Structure should be characterized by clarity and simplicity of flow 
rather than by complex branching and looping. A common type of problem is a 
program where flags are used to select overlapping sections of code to calculate 
several somewhat different functions. Modification of the code to change one 
function may unexpectedly change the other functions as well. 
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The various processes or subprograms of the model should be isolated from 
each other (modularized) so that none of the subprograms need to know the inter- 
nal structure cf the other subprograms (i.e. they should be black boxes to each 
other). With this type of structure, changes in how a subprogram does something 
will remain isolated within that subprogram and will not "ripple” throughout the 
rest of the program. 

Data flow within the program should be highly restricted so that each 
subroutine has access to Just those variables required for its operations. When 
unneeded variables are passed to a subroutine, changes in that subroutine may 
unexpectedly affect those variables, causing a problem to appear in some other 
part of the program. It can be very time-consuming to trace those problems back 
to their origins. 

Validation . Validation refers to the testing of a model on an independent 
data set. When this has been done, we have an empirical Indication as to the 
model's applicability to years and/or locations other than those for which it 
was fitted or calibrated. If a model has not been validated, then some indepen- 
dent tests should be run on the model before it is definitely considered as a 
candidate for large area yield forecasting. None of the models reviewed in this 
report have been adequately tested over a wide range of conditions and 
locations. All of the models need some additional testing. 

Timeliness . Timeliness refers to the capability of the model to produce 
projections of end-of -season yield or estimates of potential yield as needed 
throughout the growing season. For example the model should be capable of pro- 
viding a reasonably reliable yield forecast as early as it is needed. 

Subsequent forecasts should occur as needed and when significant improvements in 
earlier forecasts are possible. 

IV. Model Reviews Format 


A. Title, Author, and Reference 

This section gives the title(s) of the published report (s) on the 
model, the author(s), and the reference(s) for the publications ). 

B. Abstract or Summary 

This section consists of the abstract or sumnary from the model 
publication. 

C. Status and Applicability 

This section Indicates if (and how long ago) model development has 
been completed. If it is known, the availability and language of 
model computer code is reported. The geograhic area and environmen- 
tal conditions where the model should be applicable are Indicated. 
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D. Model Design 


This section presents the important characteristics of the model as 
reported in the manuscript. The time period used for model develop- 
ment and the time period selected for testing (if specified) are 
presented. Data requirements, which include location, meteorological 
variables, agricultural statistics and phenological Information (if 
needed), are listed including basic derived agrometeorological para- 
meters such as soil mo Is wore and potential evapotranspiratlon. 

Finally, assumptions are usually made In model development with 
regard to technological changes, weather-technology interaction, crop 
calendars (biological or phenological time as opposed to fixed or 
calendar time) and other unique features which much be considered for 
application testing. Some assumptions are also presented In Section E. 

E. Critique 

A brief surrmary of the model is presented in Section E. Section B 
(Abstract) contains a synopsis of the paper while Section E 
(Critique) is the reviewer's assessment of the capabilities and limi- 
tations of the model. It is an attempt to highlight those features 
of the model which can be identified by the reviewer to provide 
guidance for further examination In the application testing of 
the model. 

F. Author's Comaents 


This section is provided to Include any additional oorrments by the 
model authors in the review of the model. 
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A. Title, Author and Reference 


CORNF: A dynamic growth and development model for maize 

(Zea mays L.) 

by 

M. S tapper and 3. F. Arkin 


P tv gram and model documentation No. 80-2. Texas Agric. Exp. Statical, 
Te.cas A A M University, College Station, Texas. U.S.A. (1980). 

B. Summary 

A dynamic simulation model (CORNF) was developed to simulate daily 
growth and development of maize plants grovn In a wide range of 
environments. Photoperlod-temerature-genotype Interactions are 
accounted for by introduction of a maturity rating system for maize 
genotypes. Nine maturity classes are distinguished. CORNF was tested 
at nine different latitudes across the U.S.A. Observed and simulated 
values were compared for phonological stages, dry matter production, 
grain yield, yield components, leaf appearance, leaf area and water 
use. 

C. Status and Applicability 

1. Status and Availability, me model Fortran code is listed In the 
documentation and is available upon request, me model was com- 
pleted and published In 1980. Some work may be needed to adapt 
the model Fortran code to run on another computer system. 

2. Applicability, me model is designed to simulate maize pro- 
duction anywhere maize is grown by choosing the right maturity 
class, me model should be well suited for large area yield 
forecasting work because of its relatively simple Input and 
output requirements. 

D. Model Design 

me model consists of major subroutines Interacting to calculate pho- 
nological stage, leaf area development, light Interception, photo- 
synthate partitioning, ear development and water balance components. 
Growing Degree Days are used to compute changes In the development of 
the plant and are, therefore, the controlling parameter of the model. 

1. Model Development, me model was developed from material pub- 
lished In the literature and field data collected at Temple, 

Texas In 1978 and 1979. me model was calibrated with field 
data collected In 1979 at Temple (Texas), Manhattan (Kansas) 
and Swift Current (Saskatchewan, Can.). 
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2. Model Testing. The model was tested on dat* *hich had not been 
used In no del development. Some data were collected In 1979 an d 
others were published in the literature. Nine locations In the 
U.S.A. between latitudes 31 and 47 degrees were represented 

In the final test, with a total of 31 genotype-years: Texas (8), 
Missouri (1), Kansas (1), Illinois (5), Pennsylvania (4), Iowa 
(6) and North Dakota (6). 

3. Input Data Requirements. Input data for the model Include daily 
meteorological data, planting Information and soils Information. 
The daily meteorological data Include: solar radiation, 
maximum and minimum temperatures, and rainfall. Planting 

data Include: date of planting, planting density, planting 
depth, leaf area of the first leaf and latitude of the location. 
Soils data include: Initial and potential extractable water from 
each soil layer, the thickness of each layer, and coefficient 
values for soil surface tempera tore and Stage 1 soil surface 
evaporation. The maturity classes for which the model is to be 
run must be specified. 

4. Model Output. Model output consists of daily values of phenolo- 
glcal stage, leaf stage, leaf area Index, root depth, extractable 
soil moisture within the rooting depth, evapo transpiration, ratio 
of actual to potential extractable soil moisture, daily and accu- 
mulated heat units, accumulated plant dry weight, reserve 
carbohydrates, dally Increase and accumulated grain weight, per- 

- cent grain moisture, Intercepted photosynthetlcaily active 
radiation and mean daily temperature. A summary output gives 
grain yield per hectare, above ground dry matter production per 
hectare, and number and weight of kernels per plant. 

E. Critique 


This model has the level of complexity desired In a process model for 
estimation of large area yields. The major areas of physiological 
activity and some of their Interactions are modeled. As noted by the 
authors, some potential water stress effects are not modeled. These 
effects are reduction of leaf area due to water stress and delay of 
an thesis due to water stress. Overall, the model test results were 
very good. 

P. Author’s Oomnenta 

The results obtained with CORN? confirm the value of modeling crop 
phonological development when simulating crop production. The simu- 
lated total number of leaves on a plant is the essential parameter In 
the phasic development of the plant. It is used In computations to 
determine the ocurrence of tassel Initiation, ear Initiation, an the- 
sis and maturity, as well as the size of each Individual leaf. The 
leaf number is determined as a result of photoperiod-temperature- 
genotype Interactions. These Interactions are not yet fully 
understood. A maturity rating system for maize genotypes was Intro- 
duced on simplified relationships between total leaf number and 
daylength. A relationship between total leaf number and temperature 
also exists, but could not be quantified. 
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Model validation Indicated that one genotype might have to be fitted 
to different maturity classes In different environments. A par- 
ticular genotype, however, always belongs to the same maturity class 
in a given region. This gave the following distribution of maturity 
classes (1-9) over North America: from 1 In Saskatchewan to 9 in 
southern Texas and from 6 In Nebraska to 2 in Pennsylvania. 

The Growing Degree Day (GDD) concept is another source of errors. 

It neglects the effects of (a; differences between air (screen) and 
plant temperatures (e.g. water stress), (b) differences between day 
and night temperatures , (c) extreme temperatures and (d) changing 
threshold tempe 'atures with advancing age of the plant. Two other 
effects that influence GDD are accounted for In a crude way. GDD’a 
are corrected l'or daylength, based on an erqpirical relationship 
established while calibrating the model. Early In the season, when 
the growing point is close to the soil surface, dally GDD's are 
reduced under conditions when soil temperatures are expected to be 
lower than air temperatures . 

The model simulates crop production under conditions Where no serious 
nutrient dtress exists. Water stress effects are Incorporated lr che 
following processes: photosynthesis, transpiration, leaf senescence 
(leaf area after anthesla), kernel number and root growth. These 
processes could not be made very sensitive to water stress because 
water stress is derived from the simulated soil water balance and 
rooting depth which have their shortcomings. 

Dry matter and grain production were generally simulated well. Yield 
component simulations (kernel number and weight) were less accurate. 
Calculated kernel number tended to be higher than measured, espe- 
cially In high evaporative demand climates. Lower kernel weights 
simulated under conditions of limited assimilate supply, however, 
compensated for this error. Grain yields were over-estimated when 
ample assimilates were available for an over-predicted number of 
kernels. 

The model is a first version. CORN? can be Improved When more 
detailed data become available, leading to a better understanding of 
the processes Involved. 
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A. Title, Author and References 


k MODEL FOR PREDICTING CORN YIEIDS FROM CLIMATIC DATA 
by R.W. HU1 A.M. Asce, and R.J. Hanks, 

ASAE Technical Paper No. 78-4030 
Fbr Presentation at the 1978 Sunnier Meeting 
American Society of Agricultural Engineers 


B. Abstract 

Hie single most Important factor that Influences crop yields from 
one location to another, or from one year to the next, Is moisture 
availability. A better understanding of how water Influences yields is 
essential for maximizing yields through water management practices. 

Hie objective of this study was to develop a model that estimates crop 
yield as a function of moisture availability during selected periods of 
growth. 

Data Inputs Include the amount of soil water in storage at the 
beginning of the season; available soli water storage capacity for the 
root zone; and daily values of rainfall, irrigation, maximum and mini- 
mum temperatures , and specific parameters for each cultlvar that rela- 
tes phenology to daylength and/or temperature . Yield is predicted as a 
function of the relative transpiration during each of the selected 
growth periods. 

When the program is used for scheduling irrigation, the required amount 
and timing of irrigation water for any planting date is determined by 
simulating the effects of applying supplemental water In Incremental 
amounts and times. Hie "best” resultant irrigation scheduling Is Indi- 
cated for any pre-selected yield level. 

Hie program does not eliminate the need for field trials, but It can be 
used for identifying management practices that will maximize yields 
through water management or the avoidance of dry periods. Thus, field 
research can be concentrated on problem areas with resultant savings of 
time and nrney. 
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A. Title, Author and References 

Model for Predicting Plant field as Influenced by Water Use 

by R. J. Hanks 

Agronomy Journal 66(1974): 660-665 


B. Abstract 

A model has been devised to predict plant yield, both total 
dry matter and grain, as a function of water use. The model is 
simple and Inexpensive to run on a computer to determine seasonal 
yields as Influenced by Irrigation frequency and amount, rainfall, 
and soil water storage. A good fit of predicted vs. measured dry 
matter yield of sorghum ( Sorghum vulgare L. ) In Colorado, com 
( Zea mays L. ) dry matter and grain yields In Israel, and com grain 
yields In Nebraska, with various water application treatments, was 
found. A basic assumption is that the ratio of actual to potential 
dry matter yield is directly related to the ratio of actual to 
potential transpiration. Evaporation from the soil is assumed to 
decrease with the square root of time after wetting as well as with 
the stage of growth. The shape of the relative yield-water use 
curve was found to be sensitive to the evaporation and transpira- 
tion assumptions made, but Insensitive to the relation used to 
describe the influence of soil water status on transpiration. 

C. Status and Applicability 

1. Status and Availability. This model was published In 1974. 

No information is available as to further development on 
this model since 1974. If the 1974 version of the model is 
not available, it could be coded with a relatively small 
effort because of the simplicity of the model concept. 

2. Applicability. The model estimates water stress effects on 
com growth and grain yield. The model is potentially 
applicable wherever water is the dominant factor limiting 
com growth and yield. In the U.S. com belt, other factors 
appear to limit com growth to a greater degree than water in 
many years so the model may not be well suited for this 
region. The model is Intended for use in irrigation 
scheduling. 

D. Model Desigi 

The model consists of a soil moisture budget, and a function 
to estimate the contribution, during each growth period, of 
water stress to the reduction of actual grain yield below the 
potential yield level. 

1. Model Development. The model was developed on data from 
irrigation treatments in Israel, and Colorado. 
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2. Model testing. The model was tested on data from 9 irrigation 
treatments at Mead, Nebraska in 1972. The model was also tested 
on data from California. 

3. Input Data Requirements. The Input data required by the model 
include daily weather data, crop parameters, soil parameters, and 
irrigation dates and amounts. Required daily weather data con- 
sist of rainfall, and poten' ial evapotranspiration. Plant para- 
meters are planting date, root growth parameters, growing degree 
days to reach each stage from planting to maturity, growth stage 
weighting coefficients for the yield function, and values of A 
and B for each growth stage where A relates potential transpira- 
tion to potential evapotranspiration and B relates potential soil 
evaporation to potential evapotranspiration. The soil parameters 
consi”t of Initial and available soil moisture for each soil 
layer and maximum potential rooting depth. 

Model Output. The model output consists of daily values of the 
soil moisture In each layer and actual and potential E and T and 
drainage. The model also gives a final estimate of grain and dry 
matter yield. 

E. Critique 

This model may have only limited applicability In the AgRISTARS 
large area yield estimation project because it considers only water 
stress effects on yield. Other factors which may be Important 
Include intercepted light, available nitrogen, high temperature 
stress, and frost damage. 

F. Authors Coiments 


The model has been widely tested since the date of publication. It 
predicted com yields well from irrigated and salinity tests in Utah, 
Colorado, Arizona, and California (Stewart et al., 1977). It has also 
predicted the effects of water stress on spring wheat yields in Utah. 
(Rasmussen and Hanks, 1976). Further tests on com gave good pre- 
dictions in Utah (Sorensen et al., 1980 and Wenda and Hanks, 1981). 

It also predicted the Influence of irrigation on alfalfa yields 
(Hanks and Retta, 1980). A modification has been used to success- 
fully predict range production under non-irrigated conditions, but 
water stress i3 the primary factor limiting growth (Wright and 
Hanks, 1981). 

We have an operation manual which has the FORTRAN program (Retta 
and Hanks, 1980). We have work done on sensitivity in an M.S. 
thesis not yet published, indicating the model is not sensitive to 
the number of layers provided there are more than 2 and which also 
Indicates the model is not very sensitive to the water content below 
which transpiration is taken to be less than potential transpiration. 
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This model was the basis for the modification to include the effects 
of sowing date on spring wheat yi^ld (Hanks and Puckridge, 1980). 

This modification uses the same input data bit predicts leaf area 
index and dry matter accumulation on a daily basis. It has not 
been tested except under the one Australian condition. We call 
this model HAPUC. This model (as well as PLANTGRO) is available 
from R.J. Hanks in FORTRAN and BASIC code. 

Note that this model predicts relative yield as influenced by water 
stress only. Other factors that Influence yield would change the 
potential yield (no water stress yield) and would not be accounted 
for in the model. Another factor not accounted for is upward flow of 
water from a water table (sub-irrigation) . 

If you want to account for nitrogen effects and water flow up from a 
water table, you need to go to a more sophisticated model (such as 
Watts and Hanks, 1980). This model is written in FORTRAN 
and is available from D.G. Watts (Univ. of Nebraska). 


* 
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A 


Title, Author and Reference 


INITIAL VALIDATION OP A WINTER WHEAT SIMULATION MODEL 

by S. J. Maas and G. F. Arkin 

Technical Paper No. 80-4010 1980 Sumner Meeting, American 
Society of Agricultural Ehgineers. San Antonio, Texas, 

June 15-18, 1980. 

B. ABSTRACT 

A model is described that simulates the daily growth and deve- 
lopment of wheat plants based on temperature, photoperiod, soil 
moisture and population density. Information generated by the model 
during the growing season includes the length of the vegetative, 
reproductive and grain filling phases, number of productive and 
unproductive shoots per plant, and splkelet number, grain number and 
grain weight per head. As an initial test of accuracy, the mode' 1 was 
executed for comparison with phenological, morphological, and yield 
data acquired from 10 fields sown to winter wheat during the 1978-79 
growing season. All fields were located In the central portion of 
the United States — 3 in South Dakota, 3 In Nebraska, 3 in Kansas, 
and 1 In Texas. Temperature data for the simulations were obtained 
from nearby weather stations, while precipitation and solar radiation 
were measured adjacent to the fields. Soil water conditions were 
determined for each field prior to sowing to allow initialization of 
the soil water balance portion of the model. An initial assessment 
of the accuracy of the model was made by comparing observed and simu- 
lated dates of emergence, floral initiation, anthesis, and maturity 
and values of shoots/plant, heads/plant, splkelets/head, 
grains /head, and weight/grain. 

C. Status and Applicability 

1. Status and Availability. The model Fortran code is available upon 
request. The model was completed and published in 1980. As the 
publication does not specify the tyoe of computer or the dialect of 
Fortran there may be some work needed to run the model on another 
computer system. Additional information about the model program 

is available in Maas and Arkln (1980a). 

2. Applicability. The model was developed for winter wheat and has a 
vernalization (cold) requirement which must be satisfied before 
jointing can occur; however, the model may be run for spring 
wheat by changing the input vernalization coefficients. 

D. Model Design 

The model grows a representative wheat plant by simulating emergence , 
winter kill effects, transpiration , soil moisture distribution , pheno- 
logical development , leaf and tiller initiation and growth, floret 
initiation and development, and grain filling. 
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1. Model Development. The model was Initially developed from 
results published In the literature for a number of spring and 
winter wheat varieties. The model was then adjusted and verified 
against detailed winter wheat data sets obtained at Bushland, 
Weslaco, and Temple, Texas and against data from North Platte, 
Nebraska. 

2. Model Testing. The model was tested on data from 10 fields 
collected during the 1978-79 winter wheat growing season in South 
Dakota (3), Nebraska (3), Kansas (3), and Texas (1). Additional 
sensitivity testing is reported by Larsen (1981). 

3. Input Data Requirements . Input data required include planting 
date, latitude, row and wi thin-row spacing, planting depth, and 
13 varietal parameters. The model also requires soils data 
including number of soil layers to be modeled, bare soil albedo, 

U and C for stage 1 evaporation, thickness of soil layers , and 
initial and maximum available soil moisture for each layer. 
Meteorological data required Include maximum and minimum daily 
tenperatures , dally solar radiation, rainfall, and snow depth. 
Show depth is acquired to address the winter kill problem and 
for the soil moisture budget. 

4. Model Output. The model is run on a daily basis and generates 
values for growth stage, tiller number, leaf area Index, rooting 
depth, evapotransplration . soil water balance, and snow melt. A 
surrmary output Includes the number of tillers per plant, grains 
per plant, grain weight per plant, heads per area, yield per 
unit area, and yield components per tiller. 

Critique 

This model has about the level of complexity desired for large 
area yield modeling. For winter wheat, the authors noted that 
the model consistently predicted floral Initiation (approximately 
Jointing) about 2 weeks late, however they noted that the observed 
floral initiation dates might be In error. 

The model does not inlcude a photosynthesis subroutine, possibly 
because the authors felt that photosynthesis is not usually a 
primary limiting factor for wheat yield, but is itself limited by 
such factors as water stress, temperature, and mineral supply. 

However sometimes water and temperature stress limit yield 
through reduced leaf area duration, which can be modeled through 
photosynthesis and biomass accumulation. This model may 
require some revision to be suitable for spring wheat prediction 
in North Dakota but the effort may be worthwhile. Otherwise the 
leaf and tiller growth functions and the soil moisture budget may 
be suitable for Incorporation into other models. 

Authors Comments 


No additional comnents 
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A. Title, Author, and Reference 

MODELLING DAILY DRY-MATTER PRODUCTION AND YIELD OP WINIER WHEAT* 

Tbm Hodges 


B. ABSTRACT 


Applicability of many plant growth models are limited because of the 
input data requirements . Photosynthes- and respiration equations 
»<ere developed from from meteorological data that could easily be 
obtained. The single crop parameter required was leaf area Index. 
These equations were developed for winter wheat ( Trltlcum aestlvum , 
(L.)) from measurement of carbon dioxide exchange rate with field 
synthesis rates after Jointing were attributed to sink enhancement of 
photosynthesis. Respiration was estimated as a photosynthesis- 
dependent growth component and a temperature-biomass dependant main- 
tenance component . The equations predicted dry-matter accumulation. 
Head weight and yield equations were developed using predicted dally 
growth, LAI, and meteorological variables. Improvement of yield pre- 
dicted was discussed. 

*From, Photosynthesis, Growth, and Yield of Sorghum and Winter Wheat 
as Functions of Light, Temperature , Water, and Leaf Area. Chapter 
1. Ph.D Dissertation Department of Agronomy, Kansas State 
University, Manhattan, Kansas, 1978. University Microfilms 
#78-21,869. 

A. Title, Author, and Reference 

MODELING COMPONENTS OP LEAP AREA IN WINTER WHEAT* 

Jeff Baker 


B. ABSTRACT 


A leaf area index (LAI) term is used by many growth models to esti- 
mate various quantities such as photosynthesis and evapotranspiration. 
Because manual techniques for measuring LAI are extremely tedious, 
alternative methods of estimating LAI for these models are being 
sought. A model for predicting LAI was developed using the indivi- 
dual components of LAI: tillers/plant, leaves/tiller, leaf 

area/leaf, and leaf area/plant as a function of plant population, 
growth stage, water, and temperature . Equations were developed for 
winter vtfneat ( Tfrltlcum aestlvum L. em. Thell . ) to provide a dally 
estimate of each component from easily obtainable meteorological data 
and data collected on plant growth In the field. 
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The model was tested on Independent data sets from Texas, Arizona, and 
North Dakota. The model performed best when predicting LAI for fields In 
which soil moisture became limiting, but failed to match higher values of 
LAI on irrigated fields In which soil moisture did not become limiting, 
presumably because of an Inability to adequately assess the effects of 
water on average leaf size. 

* From Modeling Tiller Production and Components of Leaf Area in Winter 
Wheat as £ 'Tec ted by Temperature, Water, and Plant Population. Chapter 
l.MS. Thesis Department of Agronomy, Kansas State University, Manhattan , 
Kansas, 1982. 

A. Title, Author, and Reference 

MODELING TILLER PRODUCTION AND SURVIVAL IN WINTER WHEAT* 

Jeff Baker 


B. ABSTRACT 


Yield models requiring a leaf area index (LAI) term for the estimation 
of various quantities such as photosynthesis and evapotranspiration are 
.hindered by the fact that methods of estimating LAI in the field are 
time consuming and costly. Because the leaves of a wheat crop grow on 
tillers, as the first step in ultimately developing an LAI model, a 
plant growth model was developed to predict tiller production and 
senescence for winter wheat ( Trltlcum aestlvum L. em. Thell) on a per 
plant basis as a function of plant population, growth stage, water, and 
temperature . The equations used in the model were developed from easily 
obtainable meteorological data and data collected in the field on five 
cultivars of winter wheat hand planted on two different planting dates. 
Using a concept developed by Friend (1965) from growth chamber 
experiments, increases in tillers/plant were modeled as following the 
Fibonacci series during the vegetative phase of growth until this rate 
was limited presumably by competition and/or limiting soil moisture. In 
order to adapt this concept to the field environment, accumulated dally 
thermal units C£Tu whereXTu * £ (TMAX + TMIN)/2) with a base tem- 
perature of 0°C were substituted for the chronological time used in 
Friend's experiments. 

*From Modeling Tiller Production and Components of Leaf Area in Winter 
Wheat as Affected by Temperature , Water, and Plant Population. Chapter 
2. Ms. Thesis Department of Agronomy, Kansas State University, 

Manhattan, Kansas, 1982. 
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C. Status and Applicability 

1. Status and Availability. The tiller and leaf growth portion of the 
model was completed at the beginning of 1982. The current model card 
deck and test data set are now available but the model computer code is 
In the process of reprograi.ming for greater computer efficiency and ease 
of comprehension . The new program will probably be available sometime 
in 1982. The model is programmed In Fbrtran and could probably be run 
on another computer system with little effort. A modified version of 
the model uses satellite estimates of leaf area Index Instead of modeled 
values (Brakke and Kanemasu, 1979; Mohiuddln and Kanemasu, 1982). 

2. Applicability. The model was developed for winter wheat and will 
require son*.' adjustment, primarily of the phenology submodel to be used 
for spring wheat. The model may be suitable for use In those areas 
where solar radiation, water stress, and high tenperatures are the most 
Important factors limiting growth and yield. 

D. Model Deslgi 

The mode] grows a representative wheat plant by simulating emergence, 
photosynthesis, respiration , transpiration, soil moisture distribution, 
phenologlcal development, leaf and tiller Initiation and growth, floret 
initiation and development, and grain filling. The model is designed to 
accept measured values of soil moisture, growth stage, and leaf area 
Index throughout the season and adjust its simulation accordingly. 

1. The model has three primary parts. The first part is a set of light 
Interception, photosynthesis, and respiration equations derived from 
field chamber CO2 exchange rate measurements made on various winter 
wheat cultivars at Manhattan, Kansas from 1973 to 1980. The second part 
is a set of tiller and leaf growth functions which provides daily LAI 
values to the first part of the model. The third part Is a yield func- 
tion which accumulates dry matter production over growth periods to 
estimate kernel number, kernel weight, and grain yield. The tiller and 
leaf growth functions were developed from two plantings each of five 
cultivars of winter wheat during the 1979-80 season at Manhattan, 

Kansas. The yield function was developed from eleven plantings of four 
winter wheat cultivars during the 1976-77 season at Manhattan, Kansas. 

2. Model Testing. The first part of the model was tested for accuracy 
of biomass accumulation on winter wheat data from fourteen commerical 
fields In Riley, Finney, Ellsworth, and Colby counties In Kansas during 
1975 through 1977. The second and third parts of the model were tested 
on 24 winter wheat fields during the 1979-80 season at Bushland and 
Vernon, Texas, on 12 irrigated spring wheat fields during the 1977-78 
season at Phoenix, Arizona, and on 2 spring wheat fields during 1979 at 
Mandan, North Dakota. Jhe authors concluded that the model worked best 
on fields where water stress limited leaf growth, probably because some 
degree of water stress was present on all of the developmental data. 
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3. Input Data Requirements . The model requires Initial values 
of soil moisture for 5 layers of the soil profile as well as soil 
moisture content at field capacity and wilting point for each 
layer and two soil surface values controlling stage 1 
evaporation. The model also requires daily values of solar 
radiation , maximum and minimum temperatures, and precipitation. 

4. Model Output. The model 1'j run on a daily basis and estima- 
tes dally values of soil moisture at each level, crop biomass per 
unit area, gross photosynthesis, daytime and nighttime 
respiration, and degree of water and temperature stress. The 
model also generates a final estimate of yield and yield com- 
ponents per unit area. 

E. Critique 

This model has about the level of complexity desired for large 
area yield modeling. For use In spring wheat yield forecasting, 
the phenology submodel would have to be replaced. In the latest 
version of the model, leaf growth, tiller growth, and head growth 
functions have been included. The requirement for initial soil 
moisture content can be avoided by running the model for the 
season previous to the season When the first yield estimates are 
required. The soil related constant values have to be estimated 
for large areas and may prove to be a problem. 

F. Authors Corrments 


No additional comnents. 



A. Title, Author, and Reference 

PREDICTION OP THE INFLUENCE OP WATER, SOWING DATE AND PLANTING 
DENSITY ON DRY MATTER PRODUCTION OP WHEAT 

by R.J, Hanks and D.W. Puckridge, lAust. J. Agric. Res., 1980, 31, 

1 - 11 . 

B. Abstract 


A water balance was used to calculate dry matter yields for wheat. 

The prediction used Initial soil water, irrigation, rainfall and 
pan evaporation as inputs. Leaf area index (LAI) was estimated 
by an empirical equation and changes in LAI were determined by 
the ratio of predicted to potential transpiration and relative den- 
sity. Time of sowing influenced time of maximum LAI. Dry matter 
production was calculated from equations relating LAI and 
photosynthesis. 

The model was tested with data from wheat crops in South Australia 
which had been grown with large differences in water supply, planting 
density and sowing date between seasons. There was good agreement 
between predicted and measured production. 

C. Status and Applicability 

1. Status and Availability. The model was completed and published in 
1980. Availability of the computer code is unknovn and the paper 
does not specify the language or machine where the model was 
implemented. However, because of the simplicity of the model, it 
should not be too difficult to code again if it is not available. 

2. Applicability. The model should be applicable in any area and 
year where moisture is the primary factor limiting spring wheat 
growth and yield. It was developed for spring wheat in South 
Australia and is based on earlier models (Hanks and Rasmussen, 
Rasmussen and Kanemasu) developed for spring wheat in Utah and 
winter wheat in Kansas. 

D. Model Design 

The model is designed to estimate total dry matter production or 
biomass of a wheat crop as a function of leaf area index (LAI) 
and degree of water stress. The model estimates daily vulues of 
LAI from assumptions about the date and value of maximum LAI and 
from calculated daily transpiration. Daily photosynthesis is 
then estimated from potential evaporation, LAI, and the ratio of 
calculated actual transpiration to calculated potential 
transpiration. Dally respiration is estimated from accumulated 
dry matter. 

1. Model development. The authors used one irrigation treatment at 
the Waite Institute in Australia in 1969 as a standard year for 
model calibration. 
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2. Model testing. The model was tested on 5 Irrigation 
treatments at Waite In 1969, 5 irrigation treatments at 
Palmer, Australia In 1977, and 4 planting date treatments at 
Waite In 1970 with generally good results. 

3. Input Data Requirements. Required Input data Include planting 
date, dally or weekly potential evapotranspiration, daily preci- 
pitation or irrigation amounts, initial soil moisture, available 
soil moisture capacity with depth, and maximum LAI. 

4. Model Output. The model estimates daily Increase or decrease 
In LAI as a function of (1) days from planting, (2) planting 
density (3) date of beginning of appreciable growth (4) date 
and value of maximum LAI and (5) daily predicted and poten- 
tial transpiration . 

The model estimates daily photosynthesis and respiration from 
evaporation, predicted transpiration, predicted LAI, and pre- 
dicted dry matter accumulation. 

E. Critique 


The model should give a reasonable estimate of yield for regions 
where it is calibrated, when water stress is the major factor 
limiting yield. When factors other than water stress are limiting 
yield, the model probably shouldn’t do as well, however it may be 
possible to Incorporate it into a yield estimation system for use 
when weather data Indicates water stress is present. It may also be 
possible to modify the model by Including In it more sophisticated 
leaf growth and phenology algorithms. 

F. Author’s Conments : 

See comments unaer "A Model for Predicting Crop Yields from Climatic 
Data" by R.W. Hill and R.J. Hanks. 
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A. Title, Author, and Reference 


SOYMOD/OARDC- 

A DYNAMIC SIMULATOR OP SOYBEAN OROWIH, DEVELOPMENT, AND SEED YIELD: 

1. Theory, Structure, and Validation, by Q.E. Meyer, R.B. Curry, 
J.G. Streeter, H.J. Mederski 

Research Bulletin 1113, December 1979, 

OHIO AGRICULTURAL RESEARCH AND DEVELOPMENT CENTER 
Wooster, Ohio 

U.S. 250 and Ohio 83 South 


B. Sumnary 

SOYMQD/OARDC Is a detailed computer sinulator of the soybean plant. 
It attenpta to provide procecs description and visual description of 
the plant. An important key in this system is the breakdown and 
description of dry matter as the sum of four major entitles: struc- 
tural carbohydrate, available carbohydrate, starch, and protein. The 
mass balance system which encompasses the soli and aerial environment 
assumes that no mass is created or destroyed, but is transferred to 
or from the environeimt by the plant r.ystem. Material within the 
plant is partitioned among the morphological parts: leaf blades, 

stem-petioles, fruit, and roots. Available carbohydrate Is 
transported In the plant, and the mechanism for this has been 
described using a system of coupled partial differential equations. 

The rates of material loss from one entity subsystem nust be balanced 
by the rate of gain in another entity subsystem or returned to the 
environment. 

Simulation of soybeans or plants in general requires an adequate 
description of the internal control system. This feature of the 
plant living system has eluded the crop modeling cornnunlty for some 
time. A modest but satisfactory amount of dynamic control is pro- 
vided by a plant carbon-nitrogen balance, By expressing the role of 
carbon in the nitrogen balance and concurrently, the role of nitrogen 
in the carbon balance, the two systems are linked and work together 
as a function of rate parameters and environmental conditions. 

Carbon-nitrogen or dry natter-nitrogen ratios have little meaning to 
the total system unless they can be related to a general purpose of 
the system or specific sybsystems. This purpose is assumed to be 
mass transport and mass conversion within the plant system, in 
response to specific cues to resolve internal deficiencies at given 
locations . 
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A complete Internal control network has not been formulated. Before the 
Internal control network can be expanded, some agreement mus • be reached 
on what additional components should be modeled or related to the rest of 
the system. The system should encompass all of the components possible. 

SOYM0D/OARDC is an attempt to describe the soybean plant on the basis of 
carbon and nitrogen. This living system obviously depends on other 
nutrient components as well: phosphorous, potassium, and iron, for 

example. Future simulation efforts will address these. 

The spectrum of simulated results from this model means that simpler 
soybean models can be questioned. Statistically based crop models with 
claims of accuracy and great utility should be Judged in perspective. 
Inferential statistics were developed to aid in the testing of theories, 
but were never Intended to designate What the theory should be. 

The soybean model described in this manuscript is the original OARDC ver- 
sion to Which reference should be made. Since May 1978, versions of 
SOYMOD/OARDC have been run on computers located at Wooster, Ohio, and 
Lincoln, Nebraska, aider separate research programs. Over the past year, 
additional test and simulation runs have been performed on SOYMOD/OARDC. 
Not all of the results of these simulations are described here, but will 
be presented in future publications. 

Copies of the model will be made available to state or federal agri- 
cultural researchers provided full credit is made public. However, the 
authors assume no liability for results generated on machines outside 
their domain. SOYMOD/OARDC and more recent versions are not available as 
an extension tool, since this simulator is most suitable for research and 
teaching involving physiological processes of the soybean. 

C. Status and Applicability 

1. Status and Availability. Versions of the modr*l have been 
published in 1979 and 1980. The model is reported to be 
available to state and federal researchers. Some work might be 
needed to adapt the model to another computer system. Additional 
information on the model is available in Meyer and Curry (1981) 
and Meyer et al. (1981). 

2. Applicability. This model is more complex than is needed for the 
AgRISTARS large area yield estimation program. However, some of 
the submodels in the model may be suitably simplified and incor- 
porated into a new large area soybean yield process model. 

D. Model Design, 

1. Model Development. The major components of the model are 

described in the sum ary. The various equations in the model are 
based on an extensive study of the soybean physiology literature, 
and on numerous growth chamber studies. 
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2. Model Testing, The various model subroutines were tested againt 
field data collected in 1974, 1976, arid 1977 for Beeson soybean 
(maturity group II) grovn In Ohio. Although the test results 
were quite good, some areas were identified as needing further 
improvement Including root growth and development, nitrogen 
partitioning, Introduction of management practices, and the phe- 
nology submodel. 

3. Input Data Requirements. The input data required Include plant 
and planting data, climate data, and soils data. Plant data con- 
sist of varietal coefficients for the phenology submodel, rooting 
density characteristics for a given soil type, emergence date, 
row v,*idth, and plant spacing. Daily climate data include 
daylength, maximum and minimum air temperatures, (tew point 
tenperature, solar radiation, rainfall, and wind run. Soils data 
Include soil type, soil water retention curve, bulk density, 
hydraulic conductivity, and Initial soil water content. 

4. Model Output. Daily model output includes a complete description 
of the simulated plant with number and size of various types of 
organs and distribution of carbohydrate and protein among the 
organs. A sunmary ovtput Includes sunmarles of all meteorologi- 
cal input data, dates key growth stages are reached 'vnd a final 
description of the plant Including node number, dr^ weight 
distribution, seed weight, seed number, and total seed yield, 

E. Critique 

This model is reviewed primarily as an example of the most advanced 
crop growth and yield models. The model is probably far more sophisti- 
cated than is needed for large area yield modeling. Also some of the 
input data are not available over large areas. However, it is 
possible that some of the procedures in the model could be simplified 
for inclusion in a soybean process model. 

P. Author's Conments : 

On an IBM 370/model 3033, the model requires approximately 15 CPU 
seconds at a cost of $2.20 per run (AQNET non-prime rate). The model 
performs best on this type of computer, although thought has been 
given to a condensed second-level microcomputer version. 
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A. Title, Author, and Reference 

A Model for Predicting Soybean Yields from Climatic Data 
by R. W. Hill, D. R. Johnson, and K. H. Ryan 
Agronomy Journal 71(1979): 251-256 


B. Abstract 


The single most important factor that influences soybean ( Glycine max 
(L.) Merr.) yields from one location to another, or from one year to 
the next, Is moisture availability. A better understanding of how 
water Influences yield is essential for maximizing yields through 
water management practices. The objective of this study was to deve- 
lop a model that determines soybean yield as a function of moisture 
availability during four periods of growth. 

Data Inputs Include the amount of soil water in storage at the 
beginning of the season, available soil water storage capacity for 
the root zone, and dally values of rainfall, irrigation, maximum and 
minimum temperatures , and specific parameters for each cultlvar that 
relate phenology to daylength and temperature . 

The program predicts yield as a function of the relative transpira- 
tion coring each of four growth periods: emergence to beginning 

flowering; beginning flowering to beginning podflll; beginning pod- 
fill to end of flowering; end of flowering to maturity. 

When the program is used for scheduling irrigation, the required 
amount and timing of irrigation water for any planting date is deter- 
mined by simulating the effects of applying supplemental water in 
incremental amounts and times. The "best" resultant irrigation sche- 
duling is indicated for any pre-selected yield level. 

The program does not eliminate the need for field trials, but it can 
be used for identifying management practices that will maximize 
yields through water management or the avoidance of dry periods. 

Thus, field research can be concentrated on problem aro->c with 
resultant savings of thne and money. 


C . Status and Appl lcaMixty 

l . Status fer.d Availaoility. The 'node! was published in 1979. The 
noael was coded in Fortran IV and a source duck may be available 
upon request. If a source deck is not available it would be a 
fairly small task to code the model again because of its 
simplicity. 
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2. Applicability. The model should give reasonable results In those 
cases where water is the major factor limiting soybean yields. 
However, in many areas, high or low temperature stress are major 
limiting factors and In these areas a water stress model may not 
capture changes In yield. 

D. Model Deslgi 

The model coefficients were calibrated with yield data from a field 
study where soybean varieties from 3 maturity groups were planted at 
4 dates at each of 3 locations in Missouri in 1971 and at 1 location 
in 1972. 

1. Model Development. The model estimates relative plant dry matter 
accumulation from relative transpiration, i.e. the ratio of 
calculated transpiration to potential evapotranspiration. Growth 
stages are estimated from temperature and day length with the 
Major et al. (1975) soybean phenology model. Yield is estimated 
as a function of: relative transpiration during each of 5 growth 
stages, relative transpiration for the whole growing period, and 
a lodging factor due to excessive moisture during growth stages 2 
and 3. Calculated transpiration is reduced if the soil moisture 
budget indicates that the available water is less than l£ of the 
available water storage capacity. 

2. Model Testing. The model was tested against yield data from 9 
plantings of each maturity group in 1973 and one planting of each 
maturity group in 1976. the model had an R 2 of .96 or better 
against the test data for each maturity group. 

3. Input Data Requirements . The input data required by the model 
include planting date, maturity group, initially available soil 
moisture and available soil moisture capacity for the root zone, 
and daily values of rainfall, irrigation, maximum and minimum 
temperatures, and solar radiation. If the soybean cultures are 
not in one of the five maturity groups then growth stage dates 
are needed as the phenology submodel does not have coefficients 
for other maturity groups. 

4. Model Output. Model output consists of daily estimates of growth 
stage and soil moisture content and a final estimate of yield 
per unit area. 

E. Critique 

The model appears to be suitable for use in the large area yield 
forecasting program for areas where available moisture is the domi- 
nant factor limiting yields and where only soybean cultivars in 
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maturity groups I-V are grovei. In much of the U.S., high and low 
temperature or solar radiation are Important determinants of soybean 
yields. The model was designed to optimize Irrigation scheduling 
rather than to forecast large area yields. 

P. Authors Conments : 

The model while used to illustrate the possibility of yield changes 
due to supplemental irrigation, is not limited to this. The calibra- 
tion was completed with rainfed conditions while the verification 
data Included only 1 irrigated experiment from Missouri. An 
additional verification was obtained from a linesource sprinkler 
experiment at Kaysville, Utah, 1980. The model matched the bean 
yields quite well for the Intermediate to high irrigation treatments 
with an R 2 *88t overall treatment (uncorrected for the mean). The 
model predicted yields were considerably lower than the observed 
yields on the nonirrigated plot. 

The model uses "potential" yield levels which correspond to reaso- 
nably good management practices such as weed control, lnnoculation, 
etc. These should be adjusted to better match a particular region’s 
expected high yield potentials. The simplicity of the model makes it 
a relatively easy task to recalibrate with locally available data. 

The computer program as written can calculate evapotranspiratlon by 
any one of 5 different methods ranging from pan evaporation to the 
Penman combination equation, depending on data availability and pre- 
ference of the user. Additional crops (com, alfalfa and wheat) have 
been added to the model as of the present. 
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A. Title, Author, and Reference 

WHEAT GROWIH AND YIELD MODEL: FORTRAN PROGRAM 
J.T. Ritchie, Unpublished 

B. Abstract 
None available 


C. Status and Applicability 

1. Status and Availability. The model has not been published and 
will probably receive some further development before it is 
published. A working fortran program and test data set are 
available for testing or application. The model is intended for 
use cai winter and spring wheat. It is now being altered for 
winter and spring barley. 

2, Applicability. The model should be applicable anywhere that a 
wheat variety is grown provided that genetic coefficients are 
available for that variety. 

D. Model Design 

The model grows a representative wheat plant by simulating emergence, 
leaf growth and senescence, biomass, soil evaporation, transpiration, 
soil moisture distribution, phenological development, tiller number, 
grain number and grain filling. Each process is simulated for each 
day based on the previous values for all the processes. 

1. Model Development. The model was developed from published 
reports and published and unpublished data for winter wheat and 
spring wheat in North America, Europe, South Africa, and 
Australia. For each process, the form of the equations and the 
values of coefficients were based on physiological theory and 
experimental data from both controlled environment and field 
conditions. 

2. Model Testing. The model has been used or tested by several 
research groups on a variety of data sets with encouraging 
results. 


3. Input Data Requirements. The model requires Initial soil 

parameters, plant genetic parameters, latitude, planting date, 
and planting density. It also requires daily solar radiation, 
rainfall, maximum and minimum temperatures , and irrigation 
amounts If the crop was irrigated. 
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4. Model Output. Output from the model consists of dally values of 
soil moisture distribution, leaf and tiller number and dry 
weight, above ground plant dry weight, growth stage, grain 
number, and grain weight. The model also produces a summary out- 
put of kernel weight, grain number per unit area and yield per 
unit area. 

E. Critique 

The model is based on sound theoretical principles and the major 
plant processes contributing to yield variation are Included. 

Although the model has not been published or documented, it has bean 
tested Independently at several locations with encouraging results. 
Genetic coefficients are available for 44 varieties of winter wheat. 
Having a range of values of genetic coefficients available may 
simplify the problem of getting values of these coefficients for 
large area yield estimation applications. Like the other available 
physiological models this model has not been adequately tested so it 
is not possible to tell how accurately it will estimate yield under a 
variety of conditions. 

F. Authors Comnents : 

The model is available to any potential user in its present state. 

We only request that users communicate problems associated with its 
use for possible Improvements. The March, 1982 version has not been 
changed much from our version of a few months before then, so we are 
hopeful that the model will not be changed so frequently as it was 
during 1980 and 1981. The model is not accurate in estimating tiller 
and head number, but those numbers are not critical to yield estima- 
tion for this model. We anticipate model documentation in 1982. 
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A. Title, Author, and Reference 


A Law of the Minimum Spring Wheat Yield Model 
R.B. Cate and D.E. Phinney 
Unpublished 

B. Abstract 


No abstract available 

C. Status and Applicability 

1. Status and Availability. The model was developed In 1979 and has 
not been modified since then. Fortran programs are available to 
run the model with existing coefficients or to rederive the 
coefficients. 

2. Applicability. The model was developed and tested on the U.S. 
Great Plains spring wheat region . It has not been tested on any 
other areas so its range of applicability is unknown. 

D. Model Design . 

The model is based on Liebig’s Law of the Minimum or Law of Limiting 
Factors. Spring wheat yield in crop reporting districts (CRD's) is 
assumed to be limited by the lowest of several factors, each calcu- 
lated over several growth periods. The growth periods are calculated 
with the Robertson spring wheat phenology model using Feyerherm's 
planting model to get a planting date. The limiting factors are 
(available) nitrogen (based can estimated soil nitrogen, applied 
nitrogen, and estimated soil moisture) from planting to Jointing, (2) 
from Jointing to heading, and (3) from heading to maturity, and (4) 
average temperature from milk stage to maturity. 

Trend is accounted for by using yearly values of Feyerherm's 
relative yielding ability (VYA) factor for CRDs and by Increases 
In applied nitrogen. 

1. Model Development. The model coefficients were derived using a 
maximum likelihood algorithm. Developmental data consisted of 
meteorological data, and crop data at the CRD level. Crop data 
Included average yearly planting dates, yield values, applied 
nitrogen values, and VYA factor values for CRDs in North Dakota, 
South Dakota, Men tana, and Minnesota for 1955 to 1966. 
Meteorological data included mean daily temperatures , and rain- 
fall for estimated growth periods at the CRD level. 

2. Model Testing. The model was tested with same type of data and 
over the same regions on the years 1967 to 1976 with good 
results. 
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3. Input Data Requirements. The Input data required to run 
the model for CRD's where the model has bean derived Include 
daily temperatures and rainfall representing that CRD, nitrogen 
applied In the CRD, and average VYA for the varieties planted in 
the CRD In a given year. Values from one or several weather sta- 
tions may be used to represent a CRD, nitrogen sales data may be 
used to estimate applied nitrogen, relative yield ability may be 
derived from varietal test results weighted by acreage planted to 
each variety. The developmental and test data sets were prepared 
by Dr. Arlen Feyerherm at Kansas State University for a spring 
wheat yield model based on a least squares fitting approach. 

4. Model Output. Model output for a year and CRD consists of the 
value of each limiting factor and the maximum yield permitted by 
the most limiting factor. 

E. Critique 


The model is based can sound physiological principles. It performed 
well when tested on the 1967-1976 years in the northern U.S. Great 
Plains. The model has not been tested can any other regions so 
its geographic range Is unknown. Availability of applied nitrogen 
and varietal yielding Information may be a problem In some regions, 
but should be available for 1977 to the present for the U.S. Great 
Plains. The model should be considered for testing, as is, or should 
be refitted on more recent years if resources are available. 
Additional information on the model is available In Cate et. al 
(1979). 

F. Author's Comnants 


No additional comnents. 
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